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Abstract 

We present high resolution measurements of the thermal expansion coef- 
ficient and the magnetostriction along the a-axis of (UuGeC^ in magnetic 
fields up to 16 Tesla. From the pronounced anomalies of the lattice constant 
a occurring for both temperature and field induced phase transitions clear 
structural differences between the uniform, dimerized, and incommensurate 
phases are established. A precise field temperature phase diagram is derived 
and compared in detail with existing theories. Although there is a fair agree- 
ment with the calculations within the Cross Fisher theory, some significant 
and systematic deviations are present. In addition, our data yield a high 
resolution measurement of the field and temperature dependence of the spon- 
taneous strain scaling with the spin-Peierls order parameter. Both the zero 
temperature values as well as the critical behavior of the order parameter are 
nearly field independent in the dimerized phase. A spontaneous strain is also 
found in the incommensurate high field phase, which is significantly smaller 
and shows a different critical behavior than that in the low field phase. The 
analysis of the temperature dependence of the spontaneous strain yields a 
pronounced field dependence within the dimerized phase, whereas the tem- 
perature dependence of the incommensurate lattice modulation compares well 
with that of the dimerization in zero magnetic field. 

PACS: 64.70.Kb,65.70.+y,75.40Cx,75.80+q 
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I. INTRODUCTION 



Since the discovery of a spin-Peierls transition in the inorganic cuprate CuGeC>3 by 
Hase et al. |l| this unusual magnetoelastic transition occurring in quasi one-dimensional 
antiferromagnetic insulators has again attracted considerable attention. Compared to the 
well known organic spin-Peierls systems the structure of CuGeOs is rather simple. This 
fact and the possibility of growing large high quality single crystals allows for a much better 
study of the spin-Peierls phenomena in CuGeO"3 than in the organic compounds. 

Most properties expected from the well developed theory of the spin-Peierls transition, 
e.g. the Cross Fisher (CF) theory [§],f|], are observed in CuGe0 3 . For example elastic neutron 
and x-ray scattering show a doubling of the orthorhombic unit cell below the transition 
temperature Tsp — 14. 3K leading to two non-equivalent Cu sites in the magnetic chains. 
This lattice distortion leads to alternating Cu-O-Cu superexchange paths, i.e. alternating 
intrachain magnetic exchange constants. Also in agreement with the CF theory a gap in the 
magnetic excitations is observed ||, which scales with the structural order parameter, i.e 
the dimerization. Whereas these properties of the dimerized (D) phase of CuGeOa seem to 
be well represented by a model of spin-1/2 Heisenberg chains with a spin-Peierls transition, 
some significant deviations from this most simple treatment are present in the uniform (U) 
phase, i.e. for T > Tsp. Most strikingly the magnetic susceptibility in the U phase of 
CuGe0 3 disagrees with the temperature dependence calculated for one-dimensional spin- 
1/2 Heisenberg chains JT]||0] • There is evidence that the corrections which are necessary to 
explain the magnetism in the U phase also influence the spin-Peierls transition. It is argued 
for example that due to a frustration of the magnetic exchange in the quasi one-dimensional 
chains |||| T SP is strongly enhanced in CuGe0 3 . 

The influence of a magnetic field represents a further characteristic feature of the spin- 
Peierls transition, which can be directly compared to the different theoretical predictions. 
Due to the additional Zeeman term in the Hamiltonian the nonmagnetic dimerized phase is 
destabilized when applying a magnetic field. Consequently an additional phase with a finite 
susceptibility occurs at high magnetic fields, in CuGeOs for H>12 Tesla. An incommensu- 
rate lattice modulation has been predicted theoretically for this I phase and was recently 
observed by x-ray diffraction |J. However, the knowledge about this phase is still very 
limited, e.g. the spatial character of the incommensurate lattice modulation - domain walls 
or sinusoidal modulation - is still a subject of debate. 

The theoretical H-T phase diagrams which have been calculated with different treat- 
ments of spin 1/2 Heisenberg (or XY) chains differ significantly. A detailed experimental 
determination of the phase diagram in CuGeOs therefore allows for a test of these theories 
and future descriptions incorporating e.g. a frustration of the magnetic exchange. The the- 
ories yield different predictions for the positions of the three phase boundaries - U/D, U/I, 
and D/I - in reduced field and temperature scales. Moreover, the CF theory and the earlier 
theory by Bulaevskii et al. |10j predict a discontinuous first order D/I transition, whereas a 
continuous transition is obtained within the soliton picture [JTTJ . 

So far the H-T phase diagram of CuGeOs has been mainly studied via measurements 
of the magnetization. The phase boundaries roughly agree with the predictions of the CF 
theory fll2 |, though when analyzing the data in detail some deviations seem to be present. 
Up to now there is only little information on the thermodynamic and structural properties 
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in high magnetic fields, particularly with regard to the incommensurate phase. 

In this paper we will present a detailed study of the coefficient of the thermal expansion 
a along the a-axis in magnetic fields up to 16 Tesla fl^fl. As shown in Ref. WM there is 



a very large anomaly of a at the U/D transition. Therefore, a very detailed investigation 
of the H-T phase diagram is possible via measurements of a in external fields. Moreover, 
a represents a thermodynamic property similar to the specific heat. Therefore it allows to 
study the nature of the different transitions and the critical behavior. Furthermore, a is 
per definition a structural property, i.e. the temperature dependence of the relative lattice 
constant. Because of that the giant anomalies of a at Tsp measure the development of 
spontaneous strains (e) at the U/D transition. The spontaneous strains of the D phase in 
CuGeOs are surprisingly rather large (~ 10~ 5 ) and have been observed not only by high 
resolution capacitance methods but also by diffraction techniques |T5|--|TT| . 

A comparison with neutron diffraction data reveals that the spontaneous strains are 
proportional to the squared order parameter of the spin-Peierls transition. Thus, the high 
resolution measurements of a yield a precise measurement of the order parameter in the 
D-phase as a function of both temperature and magnetic field. In addition, we will show 
that reduced but still large spontaneous strains are present in the I phase, and thus the 
temperature dependence of the incommensurate lattice modulation can also be determined 
from our data. 

The thermal expansion measurements performed at fixed fields yield information on 
the phase transitions occurring as a function of temperature. In order to allow for the 
discussion of a complete H-T phase diagram, we will also present some measurements of the 
magnetostriction ]TjJ, i.e. the changes of the lattice constant as a function of a magnetic 



field at a fixed temperature. The magnetostriction shows anomalies at field driven phase 
transitions and thus it is possible to determine the D/I phase boundary which is nearly 
parallel to the temperature axis in the H-T phase diagram and therefore difficult to analyze 
by thermal expansion measurements. 

The paper is organized as follows. After a short description of the capacitance dilatometer 
we will give an overview of the experimental observations in Section |T|. The pronounced 
differences between the lattice constants of the three phases of CuGeO"3 are already visible in 
these raw data. The H-T phase diagram obtained from our measurements is presented and 
discussed in section [IV]. The following section [V] " Field dependence of the thermal expansion 
and the spin-Peierls order parameter' 1 is divided into five parts. Before discussing the 
influence of the magnetic field on the lattice constant in detail we will shortly repeat some 
results extractable from a in zero magnetic field. Then we are going to discuss the reduction 
of the anomaly of a at Tsp with increasing field by considering its relationship to the anomaly 
of the specific heat. The next part deals with the critical behavior of the order parameter, 
which differs for the U/D and U/I transitions, respectively. Hereafter the magnetic field 
dependence of the dimerization at low temperatures, which is found to be extremely weak 
in the D phase, and the pronounced discontinuous changes of the spontaneous strains at the 
field driven D/I transition are discussed. Finally in the last part of section [V] the temperature 
dependencies of the order parameter well below Tsp are analyzed. We find a remarkably 
strong field dependence within the D phase, whereas the temperature dependence of the 
incommensurate modulation compares well with that of the dimerization in zero magnetic 
field. 
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II. EXPERIMENTAL 



The single crystal of CuGeC>3 used for the present study was cut from a large crystal (80 
mm along the a-axis) grown by a floating zone technique ||19|| . The sample is of cylindrical 
shape with dimensions of about 6x5x8.3 mm 3 for the a-,b-, and c-axis, respectively. Various 
experimental investigations on samples prepared in this way have already been published. 
Some details of structural and magnetic properties of these crystals are described e.g. in 
Ref. 0. 

The measurements were carried out with a new capacitance dilatometer. It was originally 
designed to allow for measurements of the coefficient of the thermal expansion a = ~ • |^ (L: 
length of the sample) in fixed external fields only. Besides, it enables to measure the mag- 
netostriction, i.e. the field induced length changes AL ( H ^ at fixed temperatures. Both types 
of measurement can be performed during a single run, i.e. for exactly the same orientation 
of the crystal in the dilatometer. During the measurements the dilatometer is mounted into 
an evacuated stainless steal tube, which fits into the 40 mm bore of a superconducting 16 
Tesla magnet. Although the new dilatometer is based on our conventional one described in 
Ref. [^TJ some significant differences exist. The main difference is a thermal decoupling be- 
tween the sample and the plate capacitor, i.e. the temperature of the sample is changed, but 
the capacitor is thermally coupled to the liquid He bath. Therefore no thermal expansion 
of the capacitor itself occurs during a measurement and the capacitance changes mainly 
reflect the length changes of the sample ("small cell effect"). Moreover, the mass whose 
temperature has to be controlled is drastically reduced. It only consists of the sample itself 
and a small sample holder (~ 50 g Cu). This allows for a rapid and accurate control of the 
temperature using e.g. a software PID-technique. The temperature range is restricted by 
the maximum heater current to T<200K. 

We use platinum (Pt 103) and "Cernox CX-1050" temperature sensors (Lake Shore) for 
temperatures above and below 100K, respectively. The field dependence of the latter can be 
neglected, since its magnetoresistance is extremely small - the deviation is 2% at 2.5K for 
H = 16 Tesla and rapidly decreases with increasing T. The length changes of the sample are 
calculated from the capacitance changes measured by a temperature stabilized capacitance 
bridge (Andeen-Hagerling) with a resolution of 5 • 10 _7 pF. Thus, length changes of less 
than 0.01 A can principally be resolved. Due to mechanical vibrations etc. the resolution 
is limited to ~ 0.1 A in practice. The thermal expansion measurements are performed in a 
continuous mode, i.e. capacitance data are recorded while the temperature is continuously 
increased with a small and constant rate, usually 2 — 3mK/s. For calibrating the dilatometer 
we measured the well known [22| thermal expansion of aluminum samples of several lengths 
(4 to 8 mm). The calibration was checked by measuring a 5mm copper sample. The relative 
resolution, i.e. the scatter of the data, amounts to ~ 5 • 10 _8 /K and the deviation to the 
data reported in the literature [212 is less than ~ 1 • 10~ 6 /K for temperatures up to 200K. 
In the temperature range below 45K, which is considered in this paper only, the maximum 
deviation is even less than ~ 1 ■ 10" 7 / K - 

For the magnetostriction measurements temperature is held constant and length changes 
are detected while the magnetic field is swept from up to 14 (or 16) and back to Tesla 
with a rate of 2.5 — 4mT/s. For calibration we performed magnetostriction measurements 
on a 5mm silicon sample. The cell effect - due to magnetic impurities and eddy currents 
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which cause stresses on the dilatometer via the Lorentz force - amounts up to 250A, which 
corresponds to ~ 5 • 10~ 6 for a 5mm sample. The reproducibility of this cell effect and thus 
the accuracy of the absolute values of the magnetostriction AL/L is better than ~ 5 ■ 1CT 7 . 
The scatter in AL/L amounts to ~ 2 ■ 1CT 9 , which corresponds to O.lA as mentioned above. 



III. RESULTS 

A. Thermal Expansion 

Since some organic spin-Peierls compounds show a pronounced hysteresis of the magnetic 
susceptibility |^3],^4] as well as of the structure when changing the magnetic field at low 



temperatures, all measurements of a were performed in the field cooled mode. After applying 
the magnetic field well above Tgp, usually at T ~ 25K, the sample was cooled down to 4K 
(or 2K). Then the data were taken while the temperature was continuously increased with a 
rate of 2 mK/s. In order to check the presence of hysteretic behavior we have also performed 
some measurements in the zero field cooled mode, e.g. after applying magnetic fields of H 
= 11T and H = 14T at 4K and after decreasing the magnetic field from 14T to 11T at 4K. 
In all cases studied the data agree with those obtained from the field cooled measurements 
at the same magnetic fields. In CuGe0 3 the hysteresis of the lattice constants as well as 
that of other properties pp]- f2^ , P^ ] in the H T phase diagram is apparently rather weak and 



restricted to a very narrow field range, i.e. the data in Figs. |1] and |2| represent the behavior 
of the lattice constant a in thermal equilibrium except for the region H ~ 12. 5T (see below). 



As shown in an earlier publication [H there is a very pronounced decrease of a along 
the a-axis at the spin-Peierls transition in zero magnetic field. This decrease is related 
to a large negative uniaxial pressure dependence dTsp/dp a ~ — 4K/GPa on the one hand 
and to a spontaneous lengthening of the a-axis in the dimerized phase on the other hand. 
Fig. [I] gives an overview of the changes of the thermal expansion coefficient along the a-axis 
occurring as a function of magnetic field. A pronounced field dependence of a is present only 
below 20K, i.e. in the spin-Peierls phase. First of all the transition temperature Tsp reduces 



with increasing field, a result known e.g. from measurements of the susceptibility p8|j29|] . In 



addition, a change of the anomaly of a as a function of the magnetic field is apparent from 
Fig. [j]. In all cases the phase transition shows up by a pronounced decrease of a. However, 
note that in Fig. [I] the anomaly for H = 14T is distinctively smaller than that for H = 8T, 
whereas the sizes of the anomalies for H = 8T and H = OT are very similar. 

This nonlinear change of the size of the anomalies is related to the different low temper- 
ature phases occurring as a function of magnetic field. In H = 8T the anomaly still arises 
from the spin-Peierls transition between the U and D phases p6| , ^8| ,|9| , whereas the anomaly 



of a for H = 14T is due to a phase transition between the U the I phases. Although the 
anomaly for H = 14T is reduced in size, a clear structural change between the U and I 
phases of CuGeC>3 is inferred from these data. Moreover, this transition is strongly pressure 
dependent and leads to a spontaneous lengthening of the a-axis, similar to that in zero 
magnetic field. 

Fig. shows an expanded view of the temperature and magnetic field dependence of a up 
to 17.5K and for OT < H < 16T. With increasing magnetic field the anomaly of a is system- 
atically shifted to lower temperatures, reflecting the decrease of the transition temperature 
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with increasing magnetic field. This decrease amounts to about Tsp(O) — Tsp(H) ~ 2K for 
fields up to H = 11T and to ~ 4K for H = 16T, respectively. Considering the shape and 
size of the anomalies the curves can be clearly classified into three groups: (i) Up to 12 
Tesla the anomalies remain nearly unchanged, (ii) A very pronounced decrease of "Aa" as 
a function of the magnetic field is present in the rather small field range between 12T and 
13T. Besides the strong field dependence this second group of curves is characterized by 
additional anomalies occurring below 10K. As shown in Fig.|3] very strange thermal expan- 
sions indicating the presence of several phase transitions as a function of temperature are 
found in the entire field range between 12T and 13T, most pronounced for H = 12. 5T (see 
Fig. |3|). (iii) In the field range between 13T and the maximum field of 16T again only one 
transition is present. The sizes of the anomalies are much smaller than those at low fields 
and they further reduce with increasing magnetic field. Despite of the limited field range for 
the investigation of this third group a stronger magnetic field dependence of the anomaly 
size than in the first group (H < 12T) is apparent even from the raw data. 

The different groups of curves reflect the different kinds of phase transitions. Below 12 
Tesla there is a transition between the U and D and above 13 Tesla between the U and 
I phases. In the field region 12T < H < 13T the three phase boundaries U/D, U/I, and 
D/I meet in a tricritical point. At temperatures below this tricritical point two transitions 
are expected with increasing magnetic field (D/I and I/U). A sequence of transitions is also 
possible and present (see Fig. |J) as a function of temperature pending on the details of the 
D/I phase boundary, which is nearly horizontal in the H-T phase diagram, i.e. occurring at 
a nearly constant magnetic field. 

From the data presented so far it becomes apparent that all transitions between the three 
phases of CuGeOa lead to pronounced anomalies of the thermal expansion coefficient, i.e. 
each phase transition causes spontaneous strains. It should be mentioned that a dimerization 
alone, which is characterized by the development of alternating distances between nearest 
neighbors, i.e. an antiferro distortion, does not necessarily lead to spontaneous strains. 
Especially, within the CF theory pressure dependencies of Tsp and thus anomalies of the 
thermal expansion coefficients are obtained only, if one adds an anharmonic coupling between 



elastic degrees of freedom and zone boundary phonons [30,31 1. Such a coupling is of course no 



general property of the spin-Peierls transition and thus the pressure dependencies strongly 
differ for the spin-Peierls compounds pO|j3T[ . 



B. Magnetostriction 

Whereas the measurements of thermal expansion are a very sensitive probe of the U/D 
and U/I phase transitions, it is rather difficult to determine the D/I phase boundary since 
the temperature dependence of the latter is very small. Therefore we also performed mea- 
surements of the magnetostriction, i.e. measurements of the magnetic field induced changes 
of the lattice constant a at a fixed temperature, where the D/I-phase boundary is crossed 
almost perpendicular. 

In Fig. f| we show Aa(H)/a(H = 0) obtained with increasing magnetic field at several 
temperatures between 3 and 18K. At the lowest temperature there is only an extremely 
small magnetostriction in the D phase and at Hd/i ~ 12.5 Tesla a large jumplike decrease of 
the lattice constant a occurs reflecting the first order phase transition between the D and 
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the I phase. With increasing temperature the "jump" of the lattice constant at the D/I 
transition strongly decreases indicating that the phase transition gradually changes from a 
discontinuous to a continuous one. Note, however, that Hd/i and the total magnetostriction 
up to 14 Tesla, i.e. a(14T)— a(0T), remain roughly constant below 11K. Moreover, a large 
magnetostriction is found in the D phase in this temperature range. A comparison with the 
thermal expansion data in Fig. |2] shows that all magnetostriction curves obtained in this 
temperature range reflect differences of the lattice constant between the D and the I phases. 
The gradual change from a clearly first order transition to a continuous one can also be 
extracted from the hysteresis of the magnetostriction around Hd/i (see inset of Fig. At 
low temperatures a small hysteresis with a maximum value of about 0.2 Tesla (at T~ 3K) 
is determined when comparing the measurements with increasing and decreasing magnetic 
field. With increasing temperature the amount of this irreversibility systematically decreases 
and vanishes for T>11K. 

At temperatures between ~ 11K and 14. 5K, i.e. for T SP (0T) > T > T SP (12T) (see 
Fig. |2|) the field driven transitions are no longer between the D and I phases. The still 
very pronounced and continuous decrease of the lattice constant a up to a critical field now 
reflects the continuous phase transition between the low field D and the high temperature 
U phase of CuGeC>3. In contrast to the behavior at lower temperatures both the transition 
field and the total magnetostriction rapidly decrease with increasing temperature. The first 
observation arises from the decrease of Tgp as a function of the magnetic field (see Fig. [|), 
which implies a decrease of Hd/u with increasing temperature. The second observation 
reflects the temperature dependence of the spontaneous strain, i.e. the continuous increase 
of the structural difference between the D and U phases with decreasing temperature. 

At temperatures above 14. 5K no phase transition is found. However, a finite magne- 
tostriction is clearly observable also in the U phase. Note that in contrast to the findings 
at lower temperatures the lattice constant a now increases with increasing magnetic field. 
This magnetostriction in the U phase is not related to fluctuations of the spin-Peierls order 
parameter. It measures the magnetoelastic coupling in CuGeC>3, which is a precondition 
for the occurrence of a spin-Peierls transition ||. It should be mentioned that this positive 
magnetostriction is found in the entire temperature range studied, i.e. up to 80K. A detailed 
discussion of the magnetoelastic coupling in the uniform phase as extracted from the mag- 
netostriction and its anisotropy will be given in a forthcoming publication. The finite and 
slightly temperature dependent magnetostriction in the U phase implies a field dependence 
of the thermal expansion for T > Tsp. However, in agreement with the findings presented 
in Fig. [l| this field dependence is estimated to be extremely weak. As an upper limit for 
the difference cn(OT) — a(14T) we obtain 2 • 10~ 8 /K for 20K < T < 80K corresponding to a 
relative change of a smaller than one percent. However, this finite magnetostriction in the 
U phase has to be taken into account for a quantitative comparison of the results of thermal 
expansion and magnetostriction measurements at low temperatures (see below). 

IV. H-T PHASE DIAGRAM 

From the measurements of magnetostriction and thermal expansion it is possible to de- 
termine a complete and precise H-T phase diagram of CuGeO"3 since both quantities show 
pronounced anomalies at the different phase transitions. The anomalies of the magnetostric- 
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tion give the fields at the D/I and D/U phase transitions at different temperatures while 
those of a give the temperatures at the U/D and U/I transitions at different magnetic fields. 
It is apparent from the raw data presented so far that (i) the phase boundaries in CuGeC>3 
follow the characteristic course expected for the spin-Peierls transition and that (ii) our 
observations are in rough agreement with the findings from other properties, e.g. the mag- 



netization ||T2| , [28|] . We did, however, not observe any anomaly of the thermal expansion or 
magnetostriction at H = 8 Tesla, where Poirier et al. found evidence for an additional phase 
boundary from ultrasonic studies p6 1 . 

The A-like shape of the anomalies of a but also of the specific heat signal the 
presence of pronounced fluctuations in CuGeC>3, which have to be taken into account when 
determining the transition temperatures. In literature quite different ways have been used 



to define Tgp from the anomaly of the specific heat p3|j32|j . Ignoring the fluctuations and 
describing the anomaly as a broadened mean field step as in Ref. |33j one would obtain a 
transition temperature of ~ 14. 5K for our crystal. However, this description is apparently 
not appropriate for our more homogeneous sample. If one assumes a (nearly) symmetric 
shape of the anomaly on a reduced temperature scale and a linear regular part |[32|| , i.e. 
the behavior expected for critical fluctuations, the transition temperature is distinctively 
smaller. It is found close to the maximum of C p or the minimum of a, i.e. at 14.16K for our 
crystal. However, the experimentally observed anomalies of the thermodynamic properties of 
CuGeO"3 all show a pronounced asymmetry even close to "Tsp". Therefore the temperature 
at the minimum of a only yields a lower limit for the transition temperature. In this paper 
we have defined Tsp at the maximum slope of the anomaly yielding Tgp = 14.35, which 
lies in between the two extreme values mentioned above. We emphasize that, despite of 
this rather large uncertainty of the absolute value, the decrease of Tgp as a function of H 
is obtainable with very high precision (better than ±0.05K) from our data, e.g. from the 
shift of the minimum in a, since the shape of the anomalies does not change significantly 
(see below). 

In Fig. |5| the H-T phase diagram of CuGeO"3 for fields parallel to the a-axis is displayed. 
The transition fields measured via the magnetostriction are defined in a similar way as Tgp, 
i.e. at the maximum value of the field derivative. As visible in Fig. [|the boundaries obtained 
from the two experimental methods perfectly agree. For the transition fields between the D 
and I phases two values are given representing the hysteresis at this phase transition. The 
lower (higher) value corresponds to the transition field found with decreasing (increasing) 
magnetic field. The inset of Fig. [5] shows the D/I phase boundary on a smaller field scale. 
The hysteresis can be clearly resolved at least up to about 10K. This hysteresis and the 
shape of the anomalies of the magnetostriction (see Fig. |j) show that the field driven D/I 
transitions are of first order. The shapes of the additional anomalies in the thermal expansion 
(Fig. H), which show up very close to this phase boundary, also indicate a weakly first order 
D/I transition SIS db function of temperature. 

Similar to the experimental results found in the organic spin-Peierls compounds [Q the 
H-T phase diagram of CuGe03 does not support the soliton theories of the spin-Peierls 
transition. Most of the arguments given in Ref. ||34|| , which favor the traditional theories, 
e.g. the CF theory, in the case of TTF-BDT(Au) also hold for CuGe03. For example the 
soliton theory predicts a continuous D/I transition, whereas experimentally a first order 
transition is found. Moreover, the D/I transition is predicted to occur at g/ipHp,/! ~ 0.6E , 
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where E denotes the spin excitation gap at zero temperature fl35lj36|. Using the excitation 



gap derived from inelastic neutron scattering studies [37] , the D/I transition should occur 



at Hd/i — 8.5T, which is significantly smaller than the experimental value. 

The experimental data agree much better with the calculations of Cross [Q and the earlier 
treatment of Bulaevskii et al. |LIJ]. A simple comparison to these calculations is possible by 
considering firstly the decrease of Tsp at small magnetic fields and secondly the tricritical 
point, where the U/D, U/I, and I/D phase boundaries meet. Both theories yield in leading 

order a quadratic decrease of Tsp (H) with H, i.e. (Tsp(H)— Tsp(0))/Tsp(0) — — t- ( ^twot ) " 
The predicted quadratic field dependence is confirmed by our data up to about 6T. The 
prefactor amounts to r = 0.21K/T, which agrees with the data reported by Poirier et 
al. ||26|| . Taking into account the g-factor of 2.15 for fields parallel to the a-axis |38[] the 
theoretical treatments of Bulaevskii et al. and Cross yield r = 0.19K/T and r = 0.16K/T, 
respectively. The tricritical point is derived from our measurements at H c = 12.25 Tesla and 
T c = UK (see Fig. |). The reduced values = 0.77 and = 0.92T/K can directly 

be compared to theory. Bulaevskii et al. find ratios of 0.63 and 1.08 |TUL whereas Cross 
calculates 0.77 and 1.03 @] for the reduced tricritical temperature and field, respectively. 
Apparently, our data very well agree with the prediction of Cross, especially with regard to 
T c . The deviations from the calculations of Bulaevskii are somewhat larger. This is expected 
for quite general reasons, since more reliable results are obtained from Cross' Luther Peschel 
treatment of the one-dimensional spin excitations than from the Hartree Fock approach used 
by Bulaevskii et al.. 

In Fig. [| the phase boundaries calculated by Cross are compared to our experimental 
findings. As mentioned above there are clear discrepancies. However, the shape of the 
boundaries is very similar and there is an easy way to obtain a perfect agreement between 
the theoretical curve and the experimental U/D phase boundary. As also shown in Fig.|5| 
the experimental data follow a theoretical curve, when the reduced field scale is divided 
by 1.12 (dotted line). The "effective" magnetic field, which acts in CuGeOs, seems to be 
larger than that calculated. Cross himself already gives some possible reasons for this kind 
of deviation For example, using a too small zero-temperature susceptibility to derive 
the phase diagram causes an "effective" magnetic field that is also too small. Thus, the 
discrepancy between theoretical and experimental H-T phase diagrams might be related to 
the striking discrepancy between the measured and calculated susceptibility of CuGe0 3 in 
the U phase. The latter has been attributed to a next nearest neighbor exchange ||[7|] in 
the quasi-one dimensional chains in CuGeOa. Recently, first calculations of the H-T phase 
diagram within a model incorporating a next nearest neighbor exchange have been pub- 
lished p9f| . At the present stage these calculation do, however, not improve the agreement 
with the experimental data. 

Since particular properties of CuGe03 might influence the H-T phase diagram, a com- 
parison to the findings in other spin-Peierls compounds seems worthwhile ||40|| . As shown by 
Hase et al. flT2| , f4Tf the H-T phase diagrams of all spin-Peierls compounds roughly coincide 
on reduced temperature and field scales. However, some small but systematic differences 
are present. This is most clearly visible for the D/I phase boundary, which is located at 
larger (reduced) magnetic fields in the organic compounds than in CuGe0 3 . Despite of the 
rather large scatter of the data for the organic compounds the same trend is also found when 
comparing the D/U boundaries of the organic compounds to our data of CuGeOs. Thus, 
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an enhanced "effective" magnetic field, which destabilizes the D phase of CuGeOs, can also 
be inferred from a comparison to the organic spin-Peierls compounds. 

In order to compare the phase boundaries to the I phase one has to take into account 
a particular property of the organic spin-Peierls compounds. A pre-existing soft phonon, 
which stabilizes the commensurate lattice deformation ||10|,f4||, seems to be a characteristic 
feature of all these systems. In CuGeC>3 the experimentally observed field at the D/I tran- 
sition is rather low compared to the calculated one. Moreover, the transition temperatures 
at the U/I transitions are larger and less field dependent than predicted by theory. Thus, 
the H-T phase diagram yields no evidence that details of the phonon spectrum favor the 
D phase compared to the other phases. In agreement with this conclusion drawn from the 
H-T phase diagram no pre-existing soft phonon has been observed so far in CuGeO"3. 

V. FIELD DEPENDENCE OF THE THERMAL EXPANSION AND THE 
SPIN PEIERLS ORDER PARAMETER 

The presented data of the magnetostriction and the thermal expansion in magnetic fields 
do not only allow for an accurate determination of the phase diagram. Information on 
thermodynamic properties of the phase transition as well as information about the field 
and temperature dependence of the structure are also obtained from these high resolution 
measurements of the lattice constant a. Before we are going to discuss the field dependencies 
in detail we shortly remind on the basic conclusions drawn from the zero field data in the 
following section. In addition, all properties, whose field dependencies will be analyzed in 
the subsequent four chapters, are defined in this section. 



A. Zero magnetic field 

The main conclusions from the anomalies of the thermal expansion coefficient at Tgp 
in zero magnetic field have already been presented in a previous publication ||14|| . From 



general thermodynamic arguments an anomaly of a is expected for each phase transition 
with a finite pressure dependence of the transition temperature. This yields not only for 
mean field jumps, where the anomalies of a and the specific heat C p are related via the 



Ehrenfest relation, but also for phase transitions strongly affected by fluctuations [[42| as the 
spin-Peierls transition in CuGeOa |]T^,|32|. Assuming a finite pressure dependence of Tsp one 
expects the same temperature dependence of the leading singular parts of C p and a at the 
phase transition. Quantitatively both quantities are related by a scaling factor measuring the 
(positive or negative) uniaxial pressure dependence of the transition temperature. Thus, as 
long as one treats the anomalies of the thermal expansion coefficient (Act) and the specific 
heat (AC) in a similar way, one can determine the uniaxial pressure dependence of the 
transition temperature in the limit of vanishing pressure. The relationship reads 



<9T 



SP 



<9p 



= V„T SP .^, (1) 



where V m denotes the volume per mol. Experimentally it is found that in CuGeOs the 
similarity between C p and a is restricted to a very narrow temperature range |Tsp— T| < 0.2K, 
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where sample inhomogeneities may strongly alter the temperature dependence. Therefore 
it is not possible to extract the critical behavior of the specific heat in CuGeOa |32]| from 
the thermal expansion data. Nevertheless, Eq. [I] is valid for CuGe03, since the hydrostatic 



pressure dependence calculated from the thermal expansion and specific heat [|14]] perfectly 
agrees with the hydrostatic pressure dependence that is directly determined via susceptibility 
measurements at finite pressures 

To obtain a quantitative measure of the anomaly size we consider the largest deviation 
of the measured thermal expansion from the extrapolated behavior above Tsp (a extr ., solid 
line in Fig. [I]) Act = max(a(T < Tgp) — a extr .), i.e. the height of the asymmetric peaklike 
anomaly. Note that Aa does not correspond to a real jump of a in the sense of a mean field 
theory. Below Tsp the anomalous contribution to a, i.e. the difference between the measured 
and extrapolated temperature dependence 5a(T) = a(T) — a extr .(T), strongly differs from 
that of the specific heat. This is found experimentally, but also expected from e.g. CF 
theory. C(T < T SP ) is dominated by the magnetic excitations and therefore shows activated 
behavior at low temperatures p4| , |32| , |33|| , whereas Sa is closely related to the structural order 
parameter (Q), i.e. the dimerization of the lattice. By definition it measures the development 
of a spontaneous strain (e) in the D phase [|14|-|r6|j . 



To describe the temperature dependence of e and 6a and their relationship to Q one can 
start from the usual Ginzburg-Landau expansion of the free energy. Taking the expression 
given by Cross and Fisher ||, adding the elastic energy l/2ce 2 (c = elastic constant) and 
assuming a coupling between dimerization and lattice strains /iQ 2 e (ji = const.), which is 



quite common for structural phase transitions |45], the free energy reads: 



7 



-SP 



1 



H 



-SP 



Q 2 + ^Q 4 + /^Q 2 e+^ce 2 



(2) 



The constants 7, r], and v are related to the mean field thermodynamic properties below 
Tsp and the leading order of the decrease of T§p in magnetic fields 0. In thermodynamic 
equilibrium without external stress, i.e. for ^ = 0, it is easy to obtain the relationship 
between e and the order parameter from Eq. |2|. The spontaneous strain is proportional 
to the square of the order parameter ||T4| , [T6| and thus the anomalous contribution of the 
thermal expansion represents the temperature derivative of Q 2 : 



Q 2 (T) 



c /m s c 
-6 T = — 
H fi 



<fa(T')dT' 



(3) 



Tsi 



In Fig. ^ we show that in zero field the scaling between the spontaneous strain obtained 
from the thermal expansion coefficient and the intensity of a superstructure reflection from 
neutron scattering data f|6| is fulfilled in the entire temperature range below Tsp (see also 
Ref. |Uj). The intensity of the superstructure reflection is proportional to the square of the 
dimerization, i.e. of the structural order parameter. Moreover, it is apparent from Fig. || 
that the thermal expansion yields a high resolution measurement of Q 2 (T), since the scatter 
of these data is much smaller than that of the diffraction experiments. 

At this point we mention that it is also possible to determine the temperature dependence 
of the spin-Peierls order parameter in external fields from Equ. |3|. In this case one has to 
define the anomalous contribution of a with respect to the extrapolated behaviour of a 
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measured in the corresponding field. For quite general reasons a finite field dependence of a 
in the U phase, i.e. a finite magnetostriction, is expected due to the magnetoelastic coupling 
which is a precondition for the occurrence of a spin-Peierls transition. This magnetoelastic 
coupling also implies a magnetic field dependence of a cxtr .. However, in our analysis of 
Q(T, H) which we will present in the subsequent sections we will neglect this field dependence 
of ctextr. for two reasons. 

First, as shown in Fig. [I] the field dependence of a in the U phase is extremely small 
(not resolvable). Note that this does not contradict a finite magnetostriction, since the 
magnetic part of a does not scale with the magnetostriction itself but with its temperature 
dependence, which is very weak in the U phase as mentioned in section |111 BL Second, 



although a ex tr. is not known exactly, it is apparent from Fig. [I] that it is much smaller than 
the anomalous contribution 5a. Thus, any reasonable choice of the background as a function 
of temperature and field yields the same anomalous contribution 5a within our resolution. 
Therefore in the following we will assume a field independent a ex tr.- Up to now there are 
only few diffraction data of CuGeC>3 in magnetic fields available to check our assumption. 
The Inset of Fig. compares the spontaneous strain in a field of 10 Tesla to the intensity of 



a superstructure reflection in a field of 9.85 Tesla [46]. Although there are small differences 



in the fields as well as their orientations with respect to the crystal axes the agreement of 
the data is quite good. We mention that the same scaling factor was used for comparing 
the zero field and field data, respectively. 

Let us now turn to the discussion of the temperature dependence of the spin-Peierls order 
parameter in zero field. In Fig. [7|we compare the experimentally determined Q 2 (T)/Q 2 (0K) 
to a theoretical BCS temperature dependence (dotted line), which has been found for other 
spin-Peierls compounds [47j. Obviously, there are strong deviations in the entire tempera- 
ture range below T SP for CuGe0 3 . For T — > T SP this is of course due to the fluctuations 
of the order parameter. As a further consequence of these fluctuations the transition tem- 
perature is reduced compared to a hypothetical mean field transition temperature (Tgp). 
Since the mean field temperature dependence of Q depends on Tsp the fluctuations may 
also explain the deviations between the measured Q 2 (T) and the dotted line in Fig. [7| at low 
temperatures. Vice versa, assuming a BCS mean field temperature dependence at low tem- 
peratures one can estimate the decrease of the transition temperature due to fluctuations. A 
corresponding analysis is shown in Fig. [?|. Below about 10K the experimentally found Q 2 (T) 
follows a BCS temperature dependence with a Tg ~ 15.8K, i.e. (Tg - T SP )/Tlg ~ 10%. 
Although there is a large uncertainty and in addition to that no real physical meaning of 
Tgp, we emphasize that a quantitative comparison to mean field theories has in principle to 
consider a decrease of Tsp which is due to fluctuations. 

Close but not too close to the phase transition the temperature dependence of the spon- 
taneous strain is well described by a power law e oc (T SP — T) 2 ^ E (Fig. [7|, see also Ref. 
From a fit to our data we find 2/3 e = 0.61(5). The large error in the exponent is a direct 
consequence of the uncertainty in the absolute value of Tsp discussed above. Renormal- 
ization group theory predicts an exponent (3 = 0.325 for the order parameter of the three 
dimensional structural transition in CuGe03 (universality class 3d Ising). This prediction is 
confirmed by our data when assuming the scaling between e and Q 2 (see Eq. |3|). However, 
note that the large error bar of 2f3 e prevents a further discrimination between different three 
dimensional universality classes (XY, Heisenberg). Moreover, deviations from the power law 
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behavior of the spontaneous strain occur very close to Tgp. This is only partially due to 
sample inhomogeneities, since such deviations are also expected for a more general reason. 
As mentioned above the critical behavior of the thermal expansion coefficient for temper- 
atures very close to T SP is expected to scale with that of the specific heat. Therefore the 
critical behavior of e = J 5a should change very close to Tsp and finite values of / 5a above 
Tsp are expected. 

In the following four sections we will investigate the field dependence of the quantities 
defined above, i.e. the field dependence of: 

(i) the anomaly size Aa, 

(ii) the critical exponent 2/? £ , 

(iii) the spontaneous strain e, and 

(iv) the temperature dependence of the order parameter Q. 



B. Field dependence of Aa 

As can already be inferred from Figs. [I] and the size of the anomalies of a markedly 
decreases as a function of H. In Fig. [| Aa is plotted as a function of the magnetic field. 
The correlation between this field dependence of Aa and the phase diagram is apparent. A 
dramatic change of Aa occurs exactly at the field, where the character of the phase transition 
changes. Therefore a clear discrimination between U/D and U/I transitions is possible, when 
considering the temperature dependence of the lattice constants. Furthermore, our data 
indicate significant structural differences between the D and I phases. The much smaller 
anomalies of a at the U/I transitions mean that the spontaneous strain in the I phase is 
strongly reduced compared to that in the D phase (see below). Besides this pronounced 
change of Aa(H) at the tricritical point the data in Fig. |8] also show a smaller decrease of 
|Aa| with increasing field for both U/D and U/I transitions. Although we could follow the 
latter only within a restricted field range, our data suggest a rather strong field dependence 
of Aa at the U/I phase boundary, which amounts to about 5%/Tesla. This rather strong 
decrease, which occurs at a roughly constant T SP , indicates the presence of continuous 
structural changes within the I phase as a function of H. 

At the U/D phase boundary we find a much weaker field dependence of Aa, though 
the transition temperature decreases stronger. Up to H = 12.15 Tesla the reduction of Aa 
amounts only to about 10%. Note that the decrease of |Aa| obtained in the entire field 
range of U/D transitions is much smaller than the decrease of the specific heat anomaly 



found for a rather small field of 6 Tesla [32 . Moreover, as visible in Fig. H Aa shows a 



nonlinear decrease with increasing field similar to that of the transition temperature Tgp(H). 
Indeed the decreases of both quantities can empirically be related to each other. The ratio 



Aa/wT S p(H) (right part of Fig. H) is a constant value for all magnetic fields, i.e. Aa(H) oc 



Tsp(H). 

The size of the anomaly of a is related to the (uniaxial) pressure dependence of Tgp 
(Eq. |l|.), which strongly differs for different spin-Peierls compounds 0. Thus, Aa itself 
is not predicted by theory. However, the field dependence of Aa can be compared - at 
least qualitatively - to theoretical predictions. The universal H-T phase diagram of spin- 
Peierls compounds implies a similar universality for a single compound when studied at 
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finite pressure. This is valid at least in the limit of vanishing pressure, which is sufficient 
to discuss our data. A pressure induced change of T SP (H = OT) also affects the scaling of 
the field axis in the phase diagram. To take this into account <9Tgp/<9p has to change as a 
function of H: 



<9T SP (H,p) 



<9p 



dT SP (Q,p) 



dp 

dT SP (0,p) 



dp 



p^O 



p^O 



' Tsp(H) _ H d(T SP (H)) \ 
, Tsp(0) T SP (0) ' OH ) 



(4) 



) +T ^m) {io ^^ ] (5) 



Note that the expression in the brackets of Eq. |5], which is derived from the quadratic 
decrease of Tgp(H) (see section fV\j , is larger than 1. This means that the absolute value of 



the pressure dependence of Tgp increases with increasing magnetic field. Combining Eq. [I] 
and [| the field dependence of Aa is given by: 

Aa(H) AC(H) ( H d(T SP (H)) \ 

Aa(0) AC(0) V T SP (H) dH ) { ' 

It is apparent from this equation that one expects a rather small change of Aa as a function 
of the magnetic field. Especially, the decrease of the relative anomaly size Act (H) /Act (0) is 
expected to be much weaker than that of the specific heat. However, a detailed knowledge 
of the field dependence of the specific heat is necessary to analyze the simple scaling we find 
between v / Tsp an d Aa (Fig. ||). Simple models for AC(H), e.g. based on the free energy 
given in ||, are not sufficient. Assuming a mean field behavior AC(H) oc Tsp(H) one even 
obtains an increase of |Aa| with increasing field, which is in disagreement with the data. 
Studies of the specific heat in high magnetic fields are in progress and a detailed comparison 
of AC(H) and Aa(H) obtained at the same crystal will be given in a forthcoming publication. 
At present we conclude that the rather small field dependence of Aa for fields below 12 
Tesla shown in Fig. [| is in qualitative agreement with the field dependence estimated from 
thermodynamic relations. 



C. Critical behavior of the spontaneous strain in external fields 

In zero magnetic field the spontaneous strain e close to T$p follows a power law expected 
for the squared order parameter (see Fig. [7]). In Fig. |3| we show e as a function of the 
reduced temperature 1 — T/Tsp(H) for several magnetic fields in a double logarithmic scale. 
From these plots a power law behavior in a rather field independent temperature range 
0.7 < T/T SP < 0.95 is found for all fields studied except those very close to the tricritical 
point, where two transitions occur (see Fig. |||). Moreover, it is apparent from Fig. that 
for low magnetic fields the slope in the double logarithmic scale is roughly constant, i.e. the 
exponent 2/3 e is constant. There is, however, a pronounced difference between the low field 
data and those at 15 and 16 Tesla, where a smaller exponent is present. 

The field dependence of 2/5 e is plotted in the lower part of Fig. || As mentioned above 
and displayed in the figure there is a rather large systematic error (±0.05) of the absolute 
values due to the uncertainty in the absolute value of T$p at a given field. This uncertainty 
does, however, not affect the field dependence of 2/3 e , since the field induced decrease of 
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Tsp is extractable with much higher accuracy from our data. In other words, an alternative 
definition of Tsp in zero magnetic field causes a shift of 2(3 e to higher or lower values, but the 
field dependence of 2/3 e shown in the lower part of Fig. |^ remains nearly unchanged. For the 
U/D transitions at low fields a large exponent of about 0.6 is found, which decreases only 
slightly with increasing field. Within the error this value is consistent with that expected 
for a 3d Ising transition {(3 = 0.325) as discussed for the zero field data. 

The critical exponents of the U/I transitions are significantly smaller than those of the 
U/D transitions and show a rather strong increase with increasing magnetic field. The 
interpretation of the findings for the U/I transition is more complicated for two reasons. 
First the nature of the I phase is still discussed controversially |J and there are to our 
knowledge no theoretical predictions for the critical exponents. Second our studies are 
restricted to a field range rather close to the tricritical point where theory predicts a smaller 
exponent (3 = 0.25. This proximity might be the reason for the systematic increase of 2(5 e 
with increasing field and thus, one may argue that the small exponents we observe for all U/I 
transitions are related to tricritical behavior. However, above H c the range of small 2f3 e is 
much larger than below H c (Fig. |9]). This pronounced difference indicates a critical behavior 
within the I phase, which differs from the rather usual one found at the U/D transitions. 
In principle the small values of 2/3 e might also arise from a different strain order parameter 
coupling in the I phase. However, as we will show below, the temperature dependence of 
the spontaneous strain for high fields indicates the usual linear quadratic coupling for the 
U/I transitions, too. Irrespective of whether the critical behavior at the U/I transition is 
intrinsically characterized by a smaller /3 t or the tricritical behavior is found in a larger field 
range, the critical exponents indicate a qualitative difference between the D and I phases. 



D. Magnetic field dependence of the spontaneous strain 

From our data of the thermal expansion coefficient in magnetic fields we can follow the 
spontaneous strain as a function of magnetic field as well as a function of temperature. 
Fig. [U] presents the development of the spontaneous strain as a function of H and T. In ad- 
dition, we show the spontaneous strain as a function of temperature for some representative 



fields in Fig. |IT| The magnetic field dependence of e is dominated by the field induced phase 
transitions. At low temperatures there is only a moderate decrease of e with increasing field 
within the D phase, i.e. for H<12T. At the discontinuous D/I transition a strong decrease 
of e occurs in a narrow field range followed by a weaker decrease within the I phase. Close to 
the D/I boundary the temperature dependence of e reflects the competition between these 
phases leading to multiple transitions. 

The field dependence of e at fixed temperatures is closely related to the field depen- 
dence of the lattice constant, i.e. to our measurements of the magnetostriction presented 
in section [TTT| . However, small but significant differences occur, when comparing the field 
dependence of the lattice constant a quantitatively with e. The field dependence of e is 
always smaller than the magnetostriction at the same temperature (see Figs. |] and pp . 
These differences arise from the finite magnetostriction above Tsp. Whereas the sponta- 
neous strain is zero and thus field independent above Tsp, the lattice constant does depend 
on the magnetic field. From the definition of the spontaneous strain one obtains: 



15 



e(H, T) — e(0, T) = / (5a(H, T)dT - 5a(0, T'))dT = AL( f ' T) - AL( f ' Tp) (7) 

To>T$p 

Eq. [7] means that the field dependence of e at a temperature T is given by the difference of 
the magnetostriction at the same temperature and the magnetostriction at a temperature T 
above T SP . Subtracting the magnetostriction data at T = 20K from the low temperature 
magnetostriction leads to a good agreement between the spontaneous strain derived from 
the magnetostriction and thermal expansion data, respectively. That means that the lattice 
constant and the spontaneous strain follow variations of T and H in a reversible way in 
almost the entire field and temperature range. The only exception is a small region around 
the D/I phase boundary at low temperatures where hysteresis effects are present in the 
magnetostriction. 

Please note that the need for correcting the magnetostriction data before comparing with 
e does not imply that the magnetostrictive effects present in the U phase are also relevant in 
the dimerized phase. The correction according to Eq.|7| is just a consequence of the definition 
of the spontaneous strain. There is, however, some experimental evidence for an unexpected 
magnetoelastic coupling in the D phase. A closer inspection of the magnetostriction data in 
Fig. ^ below about 5K shows a small increase of a with increasing field in the D phase. Since 
the sign of this effect differs from the magnetostriction found at the field driven transitions, 
it is probably not related to the field dependence of the spontaneous strain. However, there 
seems to be a correlation to the magnetostriction in the U phase, which shows the same sign 
and order of magnitude. At present we have no explanation of this finite magnetostriction 
in the dimerized phase at low temperatures. Detailed investigations of the other lattice 
constants and also for crystals with different impurity concentrations are in progress to 
clarify this low temperature behavior. 

A detailed discussion of the temperature dependence of Q in external fields will be given 
in the following section. Here we will restrict the discussion of e(H) to the low temperature 
range, i.e. for T — > 0. Within the D phase the scaling between the spin-Peierls order 
parameter and the spontaneous strain is valid as discussed in section [V A| . Thus, it is 



apparent from Fig. |TT| that the field dependence of the order parameter at low temperatures 
is extremely weak. Extrapolating the data to T = our measurements are consistent with a 
field independent Q 2 (0). In other words, not only the wave vector of the distortion is pinned 
and thus field independent up to the critical field Hq n but also the amount of the structural 
distortion for T — > does not change as a function of H. This is in striking contrast to the 
strong decrease of the energy gap in the magnetic excitations due to the Zeeman splitting 
of the triplet state |37J. The scaling between structural deformation Q and the energy gap, 



theoretically expected and observed in CuGeC>3 for zero magnetic field 0J49], is obviously 



not present at finite fields. However, we emphasize that this is not in contrast to theoretical 
treatments of the spin-Peierls transition. For example, calculations within the XY model 



really yield a field independent order parameter at T = |50 . 

A similar simple interpretation of e(H) as a measure of the amount of the structural 
distortion is not possible at and above the D/I phase transition. The spontaneous strain for 
an individual field is a consequence of the strain order parameter coupling, too. According 
to our data e is reduced but still large for the U/I transitions. However, in contrast to the 
behavior at low fields the character of the lattice distortion now chane; function of 
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H. This has been studied with x-ray diffraction showing a field dependent wave vector q of 
the lattice distortion at and slightly above the D/I transition ||. In this field range it is 
impossible to correlate the field dependence of the spontaneous strain to a single quantity. 
Assuming a homogeneous incommensurate modulation of the lattice the decrease of the 
spontaneous strain may be a consequence of a reduced amplitude of the distortion and/or a 
consequence of a wave vector dependent strain order parameter coupling. Similarly, assum- 
ing a domain-wall picture, the decrease of e may be a consequence of a reduced amplitude 
of the distortion within the individual domains and/or the number of domain walls. With 
the present knowledge of the structural deformation in the I phase it seems neither possible 
to favor one of the two models nor to relate the e(H) to a structural characteristic of the I 
phase. However, note that the large anomalies at the D/I transition clearly show structural 
differences between the D and I phases. Moreover, as visible in Fig. [11] and also in the mag- 
netostriction data, there is no indication for a field independent e for T — > as we observe 
within the D phase. Obviously the structural parameter(s) determining the spontaneous 
strain in the I phase continuously change with increasing field. 



E. Temperature dependence of the spin Peierls order parameter in external fields 

The field dependencies of the quantities discussed in the last three subsections, Aa(T SP ), 
2/3 e , and e(T = const.), mainly reflect the field driven phase transitions. In particular, there 
are only minor changes within the D phase. A completely different picture is obtained 
when considering the temperature dependence of the order parameter. Within the D phase 
the temperature dependence of e well below Tsp systematically increases with increasing 
field (see Fig. |TT|). For fields of 15 and 16 Tesla, i.e. in the I phase, a much smaller 
slope is obtained. In order to analyse the temperature dependencies of the (squared) order 
parameter we will not consider the spontaneous strains plotted in Figs. |H] and [11]. Instead 
of that we will directly investigate their temperature derivatives 5a, where small changes of 
the temperature dependence of Q 2 show up more clearly. 

As a starting point of this discussion we investigate the field dependence of 5a in the 
low temperature range. Fig. |l2|(a) shows 5a at a fixed temperature of 5K as a function 



of magnetic field. It is apparent that this field dependence strongly differs from those 
discussed in the last sections. In particular, it strongly differs from the field dependence of 
the anomaly size Aa at Tsp. Within the D-phase \5a\ systematically increases with H. At 
the D/I boundary |<5a(5K)| jumps back to smaller values and for higher fields it is roughly 
constant. We emphasize that this strong field dependence of 5a(5K) in the D phase is not 
due to the decreasing Tsp(H). As displayed in Fig. |T2](b) at a fixed reduced temperature 
(Tgp/2) a very similar field dependence of 5a is present in the D phase. However, in this 
representation \5a\ is significantly smaller in the I phase than in the D phase. This difference 
can be traced back to the different sizes of the anomalies at T SP (see Fig. |8]). In Fig- |T2l(c) 
the different sizes are taken into account by considering relative values, i.e. 5a(Tsp/2) / Aa, 
as a function of H. The result is rather surprising. The zero field value roughly coincides 
with the values at very high magnetic fields, i.e. with those of the I phase. However, the 
pronounced field dependence of 5a within the D is still present. 

In the following we first analyze Q 2 (T) for different magnetic fields within the D phase. 
In order to avoid complications due to tricritical behavior and the additional anomalies in 
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the vicinity of the D /I phase boundary (see Fig. |3|) we restrict the discussion to fields below 
11 Tesla. From the data analysis given in the previous sections one concludes that the 
decrease of the transition temperature is the main consequence of the magnetic field in this 
field range, whereas only minor changes occur in all structural quantities. In particular, the 
zero temperature value of Q 2 as well as its critical behavior are roughly field independent. 
Thus, one might even expect a universal relationship between the order parameter and the 
reduced temperature t = (Tsp(H) — T) / Tsp(H) within the D phase, i.e. a simple scaling 
behavior Q(t,H) = Q(t,0). 

However, this universality is not present at all in the D phase. To demonstrate this 
we consider 5a ■ T SP (H), i.e. the derivative de(t)/dt oc <9Q 2 (t)/<9t (Fig. |13|). Obviously 
there are pronounced systematic differences of the temperature dependence of Q in different 
magnetic fields. Close to the transition temperature <9Q 2 (t)/<9t decreases with H, whereas 
at low temperatures the opposite field dependence is present. At t ~ 0.44 all curves meet in 
a single point, i.e. <9Q 2 (t)/<9t is field independent for this particular reduced temperature. It 
is not possible to improve the agreement between the curves for different magnetic fields in 
the whole temperature range by any normalization of 5a (see inset of Fig. [13]). For example, 
when investigating Q 2 (t) divided by its zero temperature value, i.e. the quantity shown for 
the zero field data in Fig. |7|, a systematic increase of the derivative with H is obtained at 
low temperatures. Thus, the BCS-like low temperature behavior of Q(T) present for H = 
(see Fig. [7]) does significantly change with increasing magnetic field. 

Instead of a field independent Q(t) we empirically find a - to our opinion - rather sur- 
prising result. The temperature dependence of the order parameter in different magnetic 



fields is universal on an absolute temperature scale. Fig. [L4] shows a plot of 5a/yTsp(H) 
versus (T — Tgp(H)): All the curves from up to 11 Tesla, i.e. in the entire field range 
of the dimerized phase, perfectly agree within the experimental resolution. Note that the 



factor l/yTsp(H), which is used to scale the y-axis, is necessary to take into account the 
slight decrease of the anomaly size Aa. We emphasize that the temperature axis for the 
different curves in Fig. [14| is not scaled but simply shifted by the field dependent transition 
temperature. To our knowledge there is no theoretical calculation of the temperature depen- 
dence of the spin-Peierls order parameter in magnetic fields, at least none that is consistent 
with our findings. For example from the exact solution of the XY-model one obtains indeed 
an increase of dC^/dt with H at low temperatures. However, the opposite field dependen- 
cies for small and large reduced temperatures as well as the crossing at t ~ 0.44 are not 
obtained |50|| . 

Now we turn to the temperature dependence of the spontaneous strain in the I phase, 
which compares well with that in zero magnetic field. At first sight any similarity between 
high field and zero field data is covered by the much smaller size of the anomalies at the U/I 
transition. However, the data in Fig. [12] already give a first hint on this similarity: After 
normalizing the 5a axis the temperature derivatives of the spontaneous strains at high fields 
nearly agree with that in zero field at Tsp/2 (Fig. |12|c). This holds at Tsp/2 but also in the 
entire low temperature range. On reduced scales the temperature dependencies of 5a in H 
= and H = 16T are very similar (see Fig. |T5|). Only close to Tsp slight differences are 
present, which have already been discussed in connection with the critical exponent 2f3 e (see 
Fig. |^). The scaling used to obtain agreement between H= and H = 16 T data is quite 
natural, since the first - considering the reduced temperature T SP /2 - takes into account the 



18 



smaller Tsp and the second the reduced anomaly size Aa and/or the reduced spontaneous 
strain in the I phase (see Figs. || and |ll]). From Fig. 15 one thus concludes that apart from 
the absolute values of both e(T = 0) and Tsp the spontaneous strain in the incommensurate 
phase corresponds to that found in the D phase at H = 0. This similarity is also found for 
the other fields with U/I transitions. 

To connect e in the I phase to the order parameter of the incommensurate modulation, 
we can follow the procedure applied for H = 0. A strain order parameter coupling causes 
a spontaneous strain and the temperature dependence of e measures that of the incommen- 
surate modulation (Eq. |2] and |3p. There are no neutron scattering measurements of the 
superstructure reflections. Thus, we can not prove a linear quadratic coupling, as we did 
for smaller fields in Fig. g. However, the similarity between zero and high magnetic fields 
strongly indicates the same linear quadratic strain order parameter coupling in the I phase. 
This means that the spontaneous strain yields a measure of the squared order parameter 
of the I phase, too. In particular, at low temperatures the order parameter of the I phase 
follows roughly a BCS mean field behavior (see Fig. |^). Up to now there are - to our knowl- 
edge - neither measurements nor calculations of Q(T) in the I phase of any spin-Peierls 
compound. Thus, at present we can neither compare our result to other findings nor judge 
whether the presented data give a possibility to distinguish between the different models of 
the structural distortion in the I phase, i.e. domain walls or sinusoidal modulation. 



VI. CONCLUSIONS 



Using a high resolution capacitance dilatometer we have investigated structural and ther- 
modynamic properties of the inorganic spin-Peierls compound CuGeC>3. The temperature 
and field dependence of the lattice constant a have been studied via measurements of the 
thermal expansion and the magnetostricton, respectively. Pronounced anomalies are found 
at all phase transitions present in the characteristic field temperature phase diagram of 
spin-Peierls compounds. Thus, clear structural differences between the three phases are 
established from our data. The very large anomalies we observe at all phase transitions 
allow for a precise determination of the H-T phase diagram of CuGeC>3. Moreover, our data 
show that all phase boundaries strongly depend on pressure, since the anomalies of the ther- 
mal expansion and the magnetostriction are related to the uniaxial pressure dependencies 
of the transition temperature and field, respectively. Summarizing our detailed comparison 
of the experimental H-T phase diagram to existing theories, one concludes that there are 
pronounced discrepancies to the soliton theory. On the other hand the predictions given by 
Cross as well as by Bulaevskii et al. are in rough agreement with the experimental results. 
This holds for the first order character of the D/I transition and also for a quantitative 
comparison of the transition fields at the phase boundaries. However, clear deviations from 
the latter theories are also present. The temperature at the tricritical point calculated by 
Bulaevskii et al., for instance, is significantly too small. With respect to this point the 
Cross theory agrees very well with the data. Moreover, the overall shape of the U/D phase 
boundary follows Cross' calculations although the field scales differ slightly. The "effective" 
field acting in CuGeOs is about 10% larger than predicted by theory. This enhancement 
indicates that particular magnetic properties of CuGeOs determine the exact positions of 
the phase boundaries. 
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Our high resolution measurements of the lattice constant do not only yield the phase 
boundaries. Properties of the ordered phase can also be extracted. From a simple treatment 
within Landau theory one obtains a scaling between the squared spin-Peierls order param- 
eter and the spontaneous strain, which is experimentally confirmed for the U/D transitions. 
Since the anomalous contribution of the thermal expansion corresponds to the temperature 
derivative of the spontaneous strain, our data yield a high resolution measurement of the 
field and temperature dependence of the order parameter. The temperature dependence of 
the order parameter is strongly affected by fluctuations. In zero magnetic field, for exam- 
ple, there are pronounced deviations from a mean field behavior in the entire temperature 
range. Close to T SP a power law as a function of the reduced temperature is found and 
the extracted exponent is consistent with the predictions for a 3d-Ising transition. In order 
to describe Q 2 (T) at low temperatures by a BCS-like behavior one has to assume a 10% 
reduction of T SP , which is due to fluctuations. 

For the discussion of the field dependence of the spontaneous strain one has to discrim- 
inate between the different low temperature phases. Up to 12 Tesla e(T, H) measures the 
dimerization, i.e. the order parameter of the D phase. There are three main conclusions 
from our data in this field range. At very low temperatures the order parameter is nearly 
field independent. The critical behavior at the U/D transition is also rather field inde- 
pendent for a wide field range. The third finding is rather surprising. The temperature 
dependence of Q 2 is universal on an absolute temperature scale, i.e. when plotting Q 2 ver- 
sus Tsp(H) — T. In contrast pronounced systematic differences are present after scaling the 
temperature axes, i.e. when analyzing Q 2 as a function of the reduced temperature T ^T ^||^ T . 
Thus, the similarity to a BCS-like behavior found for H = rapidly vanishes with increasing 
field. 

The interpretation of our findings involving the incommensurate phase is more compli- 
cated. In this field range the wave vector of the structural deformation changes with H 
and thus the spontaneous strain is not related to a single, i.e. field independent, structural 
parameter. The field dependence of e might be a consequence of a field dependent order 
parameter as well as a field dependent strain order parameter coupling. Nevertheless, the 
temperature dependence of e at a given field reflects that of the corresponding order param- 
eter. At a fixed temperature, e.g. for T — > 0, there is a strong decrease of e with increasing 
field. This decrease is most pronounced at the first order D/I transition itself, where e re- 
duces to about 40% of its zero field value. Within the I phase e shows a further decrease 
with H, which indicates continuous structural changes as a function of field. 

For fixed fields we analyzed the critical behavior of e close to T S p as well as its tem- 
perature dependence well below the transition. The critical exponents in the I phase are 
significantly smaller than those found in the D phase. Although we could investigate the I 
phase only in a rather small field range, the critical exponents indicate a qualitative differ- 
ence between U/D and U/I transitions. Apart from the region close to Tsp the temperature 
dependence of e in the I phase compares well with that of the squared order parameter for 
H = 0. Thus, the incommensurate lattice modulation roughly follows a BCS mean field 
behavior for T — > 0, similar to the dimerization in zero magnetic field. 
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FIGURES 



FIG. 1. Thermal expansion of CuGe03 along the a-axis in different magnetic fields. The solid 
line represents the extrapolated low temperature behavior a ex tr. of the U phase. 

FIG. 2. Thermal expansion of CuGeOs along the a-axis in different magnetic fields given in 
the figure. The curves are shifted by 5 • 10 _6 /K for clarity. 

FIG. 3. Thermal expansion of CuGe03 along the a-axis in magnetic fields close to the D/I 
phase boundary. The additional anomalies reflect temperature dependent transitions between the 
D and I phases and vice versa. 

FIG. 4. Magnetostriction of CuGeOs along the a-axis at different temperatures given in the 
figure. Upper panel: Discontinuous transitions between D and I phases. With increasing T (indi- 
cated by the arrows) the character of the transition gradually changes from rather strong to weak 
first order. 

Inset: Field derivatives 1/L-<9L/<9H at 3.4K measured with increasing (o) and decreasing (•) mag- 
netic field. The hysteresis, resolved up to ~ 11K, systematically decreases with decreasing T. 
Lower panel: Continuous transitions between D and U phases. The transition fields as well as the 
overall magnetostriction rapidly decreases with increasing T. Within the U phase (T=18K and T 
= 14K for H >8T) a positive magnetostriction is present. 

FIG. 5. Magnetic field temperature phase diagram of CuGeOs derived from thermal expansion 
(circles) and magnetostriction (triangles) data. The solid line shows the result of Cross' calculation. 
An almost perfect agreement to the experimental results for the entire U/D phase boundary is 
obtained by dividing the calculated field scale by 1.12 (dotted line, see text). 

Inset: Hysteresis at the D/I phase boundary obtained from the magnetostriction measured with 
increasing (closed triangles) and decreasing field (open triangles). 

FIG. 6. Comparison between the spontaneous strain (•) and the intensity of a superstructure 
reflection (o) found from neutron diffraction ^6|. Main part: Zero field data. Inset: Spontaneous 
strain in H = 10T and neutron scattering intensity in H = 9.85T. 

FIG. 7. Symbols: Temperature dependence of the squared order parameter in zero magnetic 
field as obtained from the thermal expansion. 

Solid line: Power law ((14.35 — T) 61 representing the order parameter fluctuations close to Tsp. 
Dotted line: BCS mean field behavior for Tsp = 14.35K. 

Dashed line: BCS mean field behavior for a hypothetical Tg^ = 15. 8K (see text). 

FIG. 8. Left panel: Anomaly size Aa = max(a — a ex tr) as a function of magnetic field. 
Right panel: Aa/\/TsP versus magnetic field (see text). 
The dashed vertical lines denote the field at the tricritical point. 
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FIG. 9. Critical behavior of e. Upper panel: e versus reduced temperature on a double loga- 
rithmic scale. The curves are shifted for clarity. 

Lower panel: Magnetic field dependence of the critical exponents 2(3 e . 

FIG. 10. Spontaneous strain (o) as a function of temperature and magnetic field. The phase 
boundaries between the U, D, and I phases are given by the closed symbols (•). 

FIG. 11. Spontaneous strain as a function of T for several magnetic fields given in the figure. 
The different field ranges represent the different kinds of low temperature phases. At H=12.5 T 
the competition between D and I phase close to the D/I phase boundary is seen. 

FIG. 12. Different representations of the temperature derivative 5a = de/dT oc dQ 2 /dt of the 
spontaneous strain versus magnetic field: 

a. ): <5a(H) at a fixed temperature T = 5K. 

b. ): 5a (H) at a fixed reduced temperature t = Tsp(H)/2 
a): -5a/Aa at t = T SP (H)/2 (see text). 

FIG. 13. Derivative 5a-T SP (H) oc ) versus reduced temperatures t for different magnetic 
fields given in the figure. Inset: Same data after an arbitrary normalization of the y-axis. 

FIG. 14. Universal temperature dependence of the D phase order parameter. Note, that the 
temperature axis is only shifted but not scaled by Tsp(H). The division of 5a by ^/Tsp(H) takes 
into account the slight decrease of the anomaly size with H (see text). 

FIG. 15. Comparison between the derivatives Sa-Tsp(H) oc versus reduced temperature 

t in H = (o; left y-scale) and H = 16T (•; right y-scale). 
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